The Nuclear Ring in the Barred Spiral Galaxy IC 4933 



Stuart D. Ryder^'^ , Samuel M. IllingworthP , Robert G. Sharp^ and 
Catherine L. Parage'" 

^ Anglo- Australian Observatory, P.O. Box 296, Epping, NSW 1710, Australia 

^ Department of Physics & Astronomy, University of Leicester, Leicester, LEI 7RH, UK 

'^Research School of Astronomy and Astrophysics, Australian National University, ACT 2611, 

Australia 

^ Email; sdr@aao.gov.au 



Abstract: 

We present infrared imaging from IRIS2 on the Anglo- Australian Telescope that shows the barred spiral 
galaxy IC 4933 has not just an inner ring encircling the bar, but also a star-forming nuclear ring 1.5 kpc 
in diameter. Imaging in the u band with CMOS on Gemini South confirms that this ring is not purely 
an artifact due to dust. Optical and near-infrared colours alone however cannot break the degeneracy 
between age, extinction, and burst duration that would allow the star formation history of the ring to 
be unraveled. Integral field spectroscopy with the GNIRS spectrograph on Gemini South shows the 
equivalent width of the Pa/3 line to peak in the north and south quadrants of the ring, indicative of a 
bipolar azimuthal age gradient around the ring. The youngest star-forming regions do not appear to 
correspond to where we expect to find the contact points between the offset dust lanes and the nuclear 
ring unless the nuclear ring is oval in shape, causing the contact points to lead the bar by more than 
90°. 

Keywords: stars: formation — galaxies: evolution — galaxies: individual (IC 4933) — 
galcLxies: nuclei — galaxies: starburst 



1 Introduction 



High spatial resolution multi-wavelength imaging has 
revealed a host of intriguing features close in to the nu- 
clear regions of barred spiral galaxies, includi ng star 



l2004l : iHeller et al] 120071 '). Thus, we see that nuclear 
rings, nested bars, and AGN are all fundamentally 
inter-related, and yet we ha ve surprisingly little data 
on their rela tive frequencies |Knapenll2005t ). or causal 
connections (|Laine et al.ll200a) . 
burst rings, dust spirals, and nested bars (e.g., Rvd er fc KnapenI 



The total sample of nuclear rings of sufficient angu- 
lar size to be resolved with ground-based seeing-limited 
spectrographs or aperture synthesis int erferometers is 
still quite limited (jComeron et al.ll2010l ). Galaxies ex- 
hibiting an inner ring encircli ng their bar very ofte n 
harbour a nuclear ring as weU (|Buta fc Combeslll996l l. 
The nearly face-on barred spiral galaxy IC 4933 is 
listed in th e Catalog of Southern Ringed Galaxies (CSRG; 
'Buta| [l995l ) as having an inner ring of 47 x 31 arcsec 
across. Table [T] summarises some other properties of 
this galaxy. In this paper, we report the discovery of 
an apparent nuclear ring in IC 4933 from near-infrared 
(NIR) and optical imaging, and provide spectroscopic 
confirmation of its nature. In Sections 12.11 and 12. 21 we 
describe the image acquisition and reduction proce- 
dures, and analyse the colour maps. We present follow- 
up NIR Integral-Field Unit (IFU) spectroscopy in Sec- 
tion l2.3l which enables us to reconstruct the recent star 
formation history. Section [3] puts this in the context of 
an emerging picture of sequential star formation being 
uncovered within nuclear rings. 



11999 ; Martini et al. 2003 ; Lainc et al. 2002) . Early nu- 
merical simul ations of dynamical evolution in barred 
galaxies (e.g.. lSchwarj[T98ll ; lAthanassoula et al.|[i982l ) 
appeared to link circumnuclear starburst rings with 
the inner Lindblad resonance (ILR) , the radius at which 
the bar pattern speed flp is equal to Q,—k/2, where Q, is 
the mean angular frequency, and k the radial epicyclic 
frequency. More recent hydrodynamic simulations by 
iRegan fc Teubenl (120031 ) have argued that star-forming 
nuclear rings actually result whenever a certain frac- 
tion of gas clouds find themselves on elongated orbits 
aligned perpendicular to, rather than parallel with, the 
major axis of the bar, and that furthermore these rings 
shrink with time. 

As atomic gas in the circumnuclear regions is gen- 
erally scarce, the bar (or some form of oval distortion) 
is thought to play a significant role in 'funneling' gas 
inwards via shocks to feed the competing needs of a 
star- forming ring, and/or an Active Galactic Nucleus 
(AGN) component. When gas is trapped near or be- 
tween resonant orbits, and is consumed in star forma- 
tion, the AGN will potentially be starved and forced to 
go dormant. One way to overcome this 'logjam' is to 
create a secondary, or 'nested' bar, which can disrupt 



the resonances of the main bar (^Englmaier fc Shlosma: 
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Table 1: Basic data for IC 4933 



Parameter 



Value 



Source 



Hubble Type 
RA (J2000) 
Dec (J2000) 
Vhoi (km s"i) 
D (Mpc) 

Scale (pc arcsec"^) 
Inclination (°) 
Extinction Ab (mag) 



SB(rs)bc 
20^03™29!04 
-54°58'47'.'8 
4915 
66 
313 
37 
0.199 



1 
2 
2 
3 

4 
4 
5 



Sources: f ll lButalll99i (21 ISkrutskie et al.l [2006t 
r3'l iMever et all 20041: (4) NED; (5) 
iLauberts fc Valentiin|| l98a f6) ISchlegel et al.l 



1998 



2 Observations and Results 
2.1 Infrared Imaging 

Nuclear rings are seriously under-represented in the 
CSRG, not only because the central regions of many 
galaxies are overexposed on the UK Schmidt Telescope 
J survey films on which the CSRG is based, but also 
due to the camouflage effect of dust within the bar re- 
gion. Imaging individual galaxies at NIK wavelengths 
has a number of advantages in this regard: 

• The extinction due to dust at 2 /i m is only one- 
tenth that at optical wavelengths (jCardelli et al.l 
119891 ). 

• The rapid evolution from blue to red super giants 
after the flrst 5-10 Myr of a starburst ( Leit herer et 
11999') moves the peak of the spectral energy dis- 
tribution into the infrared. Many of these cir- 
cumnuclear star formation events are not visible 
at all optically in this period. 

• The underlying infrared background is a much 
better tracer of the dominant stellar population 
in the bar and nuclear region, allowing a bet- 
ter estimate of the mass-to-light ratio, and ul- 
timately detailed modeling of the gravitational 
potent ial which can be comp ared with velocity 
fields (jButa fc Combesll2000l ). 

The existence of a circumnuclear ring in IC 4933 
was initially suspected from the difference of J and Ks 
images obtained with the Infrared I mager and Spectro- 
graph 2 (IRIS2: lTinnev et aljr2004l ') instrument on the 
3.9 m Anglo-Australian Telescope (AAT) in June 2003 
under marginal conditions. The imaging presented 
here was obtained on 2006 Sep 4 UT in photomet- 
ric conditions, with seeing of 1.0 arcsec or better. The 
IRIS2 detector is a 1024 x 1024 Rockwell HAWAII-1 
HgCdTe array with a pixel scale of 0.45 arcsec pixeP^, 
resulting in an instantaneous field-of-view (FoV) of 
7.7 X 7.7 arcmin. 

Although the optical extent of IC 4933 is only 
~30% of the IRIS2 FoV, we took a conservative ap- 
proach and acquired matching observations of adjacent 
blank sky to track changes in the background level and 



Right Ascension (J3000) 

Figure 1: IRIS2 J-band image of IC 4933. 
The corresponding surface brightness in units of 
mag arcsec"^ is indicated by the colour bar at 
right. The inset box outlines the region shown 
in the colour index maps in Figures [5] and [31 



illumination pattern (see IVaduvescu &: McCalil (|2004f ) 
for a discussion of the merits of this strategy). Five 
jittered observations (20 arcsec offsets) of the target 
galaxy were bracketed and interleaved with six jittered 
observations of the sky 10 arcmin south. At each ob- 
aU ject or sky jitter position, a single 30 sec integration 
was recorded (in J) or 3 x 10 s integrations were aver- 
aged (for H or Ka). This pattern was repeated twice, 
for a total on-source exposure time of 5 min in each of 
the J, H, and Ks filters. 

The data reduction was carried out using the ORAC- 
DEQ pipeline within the STARLINK package and the 
CHOP_SKY_JITTER recipe. Pre-processing of aU raw 
frames included subtraction of a matching dark frame; 
linearity and inter-quadrant crosstallQ correction; and 
bad pixel masking. The first 6 sky frames are offset 
in intensity to a common modal value, then a fiatfield 
is formed from the median value at each pixel. The 
6 sky frames, and 5 object frames are flatfielded, then 
the modal pbcel values of the two sky frames brack- 
eting each object frame are averaged and subtracted 
from that object frame to account for the sky back- 
ground. A correction for astrometric distortion inter- 
nal to IRIS2 is applied by resampling the object im- 
ages, then the actual spatial offsets between images are 
computed using point sources common to all images. 
The 5 object images are mosaiced together by applying 
offsets in intensity to the registered images to produce 
the most consistent sky value possible in the overlap 
regions. A new fiatfield and mosaic is constructed for 

^http://www. oracdr.org/ 

^ htt p : / / www .eso.org/~gfingcr /hawaii_ 1 Kxl K/ 
crosstalk_rock/crosstalk.html 



www.pu blish . csiro. au /journals/ pasa 



3 




aO^Oa^aCO" 29.5' 29.0' 2B.5" 28.0" 

Right Ascension (J3000) 




20°03"30.0= 29.5" 29.0" 28.5= 28.0" 

Right Ascension (J3000) 



Figure 2: IRIS2 (J — H) colour index map of the 
inner bar region of IC 4933. The colour range dis- 
played is indicated by the scale at right. 



Figure 3: IRIS2 {H — K) colour index map of the 
inner bar region of IC 4933. The colour range dis- 
played is indicated by the scale at right. 



the second set of 6 sky and 5 object frames, then the 
two mosaics are in turn registered and co-added to 
form a master mosaic. 

Point sources detected in the master mosaic have 
been cro ss-referenced against the 2MASS Point Source 
Catalog (|Skrutskie et al. 2006) and a new astrometric 
solution computed, while photometry in a 10 arcsec 
aperture has enabled the photometric zero-points to 
be determined in each filter, after transforming from 
the 2MASS JHKs system to the Mauna Kea Obser- 
vatories (MKO) JHK system (Rvder 2007). 

Figure [1] shows the J-band (1.25 /^m) image of 
IC 4933 on a logarithmic scale. The bar, inner ring, 
and two stubby spiral arms are evident. The bright 
compact nucleus is well-fit by a Gaussian with FWHM 
of 2.0 arcsec, but shows no obvious indication of un- 
usual circumnuclear structure. Aligning and subtract- 
ing the calibrated surface brightness images to form 
colour index maps tells a different story. Figure [2] 
shows the (J — H) colour of the central 18 x 23 arcsec 
region, as outlined in Figure [T] An almost perfectly 
circular ring 5 arcsec (1.5 kpc) in diameter stands out 
by virtue of its colour; whereas the underlying stellar 
population in the bar and bulge has (J — H) ~ 0.65, 
the ring has (J — H) > 0.7 and is as red as 0.8 in the 
northwest quadrant. The {H — K) map in FigureOhas 
less contrast, but confirms the existence of this feature. 
The western half of the ring has (H — K) > 0.43, reach- 
ing 0.48 in the northwest quadrant, while the back- 
ground colour averages 0.39. 

In Figure [4] we plot the (J — H) colour against the 
{H — K) colour averaged over a 2 x 2 pixel area at the 
4 cardinal points around the ring, as well as the mean 
colours of the background stellar population. The red- 
dening vector corresponding to 2 mag of extinction in 



the V-band is also pl otted. We have used the Star- 
burst99 models from iLeitherer et al.l (|1999| ) to com- 
pute the evolutionary track through this two colour 
diagram of a stellar population having an Initial Mass 
Function with a Salpeter-like slope {a — 2.35) between 
1 and 100 Mq, and super-solar metallicity {Z — 0.04) 
typical of the bulges in Milky Way-like systems such 
as I C 4933. To encompas s most star formation histo- 
ries, [Leithere7^r|al] (|l999l ) consider two limiting cases: 
an 'Instantaneous' burst of star formation (IB), and a 
'Continuous' constant rate of star formation (CSF). 
The evolution over 250 Myr for these two scenarios is 
marked on Figured Both start off with almost identi- 
cal colours, then evolve rapidly bluewards for the first 
5-7 Myr until reversed by the evolution of massive blue 
stars into red supergiants. The destruction of these 
red supergiants in core-collapse supernovae allows the 
colours to again evolve bluewards after ~9 Myr and 
~40 Myr for the IB and CSF models, respectively, but 
now with (J — H) being 0.5 mag redder than on the 
initial 'blue loop'. 

At no stage do these evolutionary tracks achieve 
colours as red as those measured around the ring. While 
the colours of the background stellar population are 
similar to those reached by the IB model in the brief 
period 8-20 Myr with almost no reddening, this is an 
extremely unlikely scenario as the amount of nuclear 
dust highlighted in Figure [5] indicates. Rather, being 
more like the intermediate age population from contin- 
uous star formation marked by the triangle in Figure^ 
the stellar background light in the bar and nuclear re- 
gion would experience Ay ~ 2 mag of extinction as 
indicated by the reddening vector. The ring colours 
are consistent with an IB population at just 5 Myr old 
with up to Av = 6 mag of extinction, with the re- 
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Figure 4: Two colour diagram for 4 specific foca- 
tions around the circumnuclear ring in IC 4933: 
(N)orth, (W)est, (S)outh, and (E)ast. The large 
open circle indicates the spread of colours from 
unresolved stars in the bulge and bar. The ar- 
row indicates the effect on the NIR colours from 
reddening equivalent to 2 mag of extinction in the 
T^-band. The solid line is the evolutionary track 
of a stellar population formed f rom an instan- 
taneou s burst (IB) according to iLeitherer et al.l 
( 19991 ). reaching the open square after 250 Myr. 
The dashed line ending at the open triangle is the 
evolutionary track of a stellar population undergo- 
ing continuous star formation (CSF) over 250 Myr. 



quired reddening dropping to 0.5-1 mag for an age of 
9 Myr. The implied range of extinction from the CSF 
models is just Ay = 1.5-2.5 mag, but the allowable 
range of ages is much greater at 10-100 Myr. It is this 
degeneracy between extinction, age, and even the du- 
ration of star formation that requires us to seek better 
diagnostics than broadband NIR colours (Section l2.3p . 

2.2 Optical Imaging 

As Figure[3]indicates, an alternate scenario is that the 
ring is simply dust seen in silhouette, reddening and 
obscuring the background population by the equivalent 
of between 0.3 and 1.4 mag of extinction in V . While 
linear dust lanes and nuclear dust spi rals are rel atively 
frequent in spiral galaxy nuclei (Mar tini et aL| [2003). 

and compa ct nuclear du st disks have b een seen around 

some AGN (iJaffe et al.l[l99S : van der Marel fc van den BoschI 
Il998l '). such a kiloparsec-scale coherent dust ring has 
not been seen elsewhere. 

To test for the presence of a dust ring, we need a 
much larger extinction contrast than that offered by 
NIR colours. Images of IC 4933 in the Sloan Digital 
Sky Survey (SDSS) u , g' , and r' filters were obtained 
with the Gemini Multi-Object Spectrograph (CMOS; 
iHook et al] 120041 ') attached to the Gemini South Tele- 
scope as part of programme GS-2008A-Q-204 (PI: S. 
Ryder) in the "Poor Weather" queue (seeing > 1.1 arc- 
sec in non-photometric conditions with > 2 mag of at- 
mospheric extinction) on 2008 July 7 UT. Only the 
u' data is presented here. Three exposures of 690 s 
each were obtained, with each exposure offset by 8 arc- 
sec spatially to allow filling-in of the inter-CCD gaps 
in CMOS. On-chip binning yielded a pixel scale of 
0.292 arcsec pix"'^. 

The data were reduced and combined using VI. 9.1 
of the gemini package within iRAlj^. A master bias 
frame (constructed by averaging with 3-a clipping a 
series of bias frames) was subtracted from all raw im- 
ages in lieu of overscan fitting and subtraction. Images 
of the twilight sky were used to flatfield the images, 
then the dithered galaxy images were registered and 
averaged together with the imcoadd task to eliminate 
the inter-CCD gaps, bad pixels, and cosmic rays. No 
absolute flux calibration of this data is available, but 
for the purposes of constructing colour index maps the 
final image has been converted to a logarithmic scale. 

This u image has been scaled and aligned to match 
the IRIS2 K image using the geomap and geotran tasks 
before the latter was subtracted from the former to 
produce a {u' — K) colour index map which spans a 
factor of 6 in wavelength. This map has been unsharp 
masked by subtracting off a lightly smoothed version of 
itself, to help accentuate any apparent dust lane struc- 
ture. Figure [5] shows this difference in optical versus 
NIR morphology, together with a contoured represen- 
tation of the ring seen in (J—H). One can quite clearly 
see dust lanes offset either side of the bar major axis 



^IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 
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Figure 5: Contours of {J~H) colour in the nucleus 
of IC 4933 from IRIS2, superimposed on a map of 
the {u' — K) colour index in bar and nuclear region 
from the difference of CMOS and IRIS2 images. 
An unsharp mask has been applied to the {u' — K) 
image in order to highlight the two dust lanes 
which trace the leading edge of the bar and encircle 
the nuclear ring seen in [J — H). The contours cor- 
respond to {J -H):^ (0.71, 0.73, 0.75, 0.77, 0.79), 
while the bluest [u' — K) colours appear white in 
this figure. 



and sligh tly concave towards i t. As shown in the simu- 
lations of lAthanassoulal (|l992h . such offset curved dust 
lanes mark the loci of shocks at the leading edge of a 
relatively fat bar. These dust lanes curl around the 
nucleus and encircle it, though this dust ring is notice- 
ably larger (~ 8 arcsec diameter) than the NIR ring 
found in Section [2.11 Within this dust ring sits three 
blue 'hot spots', two of which are concentric with the 
NIR ring and which in fact coincide with the bluest 
(J — H) sections of the NIR ring. 

2.3 IFU Spectroscopy 

Near-IR spectroscopy can place much tighter constraints 
on the ages and star-forming timescales of individ- 
ual c lusters (|Rvder et al.l I2OOII : lAlonso-Herrero et al.l 
l200lh . The combination of multiple, independent age 
diagnostics (including the Br7 2.16 /im equivalent width; 
the CO(2,0) 2.3 iiva spectral index; or the [Fe 11] 1.26 /im 
luminosity) can in principle discriminate the IB from 
the CSF models of iLeitherer et all (| 19991 ) . as well as 
their relative ages. Since most of these diagnostics in- 
volve the ratio of an emission/absorption feature to 
their adjacent continuum emission, they are far less 
affected by differential extinction than the broadband 
colour indices in Secti0n r2.ll 



Spectroscopy of the circumnuclear ring in IC 4933 
was conducted in queue mode with the Gemini Near- 
InfraRed Spectrograph (GN IRS: Elias et al...2006') with 
its Integral Field Unit aFU: lAmngton-Smith et al.ll2006l ') 
attached to the Gemini-South Telescope for program 
GS-2005A-Q-25 (PI: S. Ryder). The GNIRS IFU slices 
a 4.8 X 3.15 arcsec FoV into 21 slices 0.15 arcsec wide 
by 4.8 arcsec long and reformats it into a single 'long 
slit', with each of the 21 slices projecting to 32 0.15 arc- 
sec pixels on the 1024 x 1024 Aladdin InSb array. The 
31.7 lines mm^^ grating and short camera disperse the 
stacked slices at a resolving power A/AA ~ 1700, with 
separate grating settings spanning much of the J-band 
(1.06-1.45 and the /C-band (1.90-2.56 ^im). 

As the ring just overfills the IFU FoV, four adja- 
cent pointings were observed, effectively covering each 
quadrant of the ring in turn. Each pointing was ob- 
served in 4 X 250 sec integrations, with 32 reads per in- 
tegration delivering a read noise of 7 electrons. Each of 
the 4 integrations was dithered by ±0.3 arcsec (2 pix- 
els) in both axes to compensate for reduced throughput 
in the two outer and one of the central slices, and to 
provide 2 pixels of overlap between pointings. Further- 
more, each integration was interleaved with equivalent 
integrations of blank sky 1 arcmin east, to enable re- 
moval of the telluric OH and thermal emission. Obser- 
vations of the FO V stars HD 176193 and HD 200767, 
nodding between each half of the IFU, were made at 
comparable airmass to the IC 4933 observations to cor- 
rect for telluric absorption. Observations were con- 
ducted over 5 nights in March and April 2005 when 
conditions were photometric, in seeing of 0.8-1.0 arc- 
sec. Spectra of an Argon arc lamp provided the wave- 
length calibration, while observations of a 1100 K 'grey 
body' or quartz halogen lamp facilitated the correction 
for pixel-to-pixel sensitivity variation. 

The IFU data has been reduced using the gnirs 
package within iraf. Images taken in 'Bright Objects' 
read mode (i.e., telluric standards, flatfields, and arcs) 
were found to exhibit a ripple pattern noise along each 
column, having a period of 17 pixels and amplitude 
~260 electrons, well in excess of the 38 electrons read 
noise expected in this mode. We have been able to 
remove this ripple pattern by median filtering the un- 
exposed columns, then subtracting this from every col- 
umn in the image. The 32 columns corresponding to 
each of the 21 slices are extracted into individual im- 
age planes. Each of the 21 object slices has the match- 
ing nodded object or blank sky slice subtracted from 
it, and is then divided by the normalised, smoothed 
fiatfield lamp slice. Tracing of an artificial 'pinhole' 
spectrum allows each slice plane to be straightened in 
the spatial direction, while interactive arc line identifi- 
cation and fitting is used to rebin each slice on a linear 
wavelength scale. 

Subsequent processing has been carried out using 
purpose-written idl routines. For the IC 4933 obser- 
vations the 21 slices are stacked into a 3D datacube of 
32 X 21 spatial pixels x 1022 spectral elements. The 
4 datacubes from each of the dithered pointings cov- 
ering each quadrant are registered and combined into 
one datacube per quadrant. The 4 quadrant cubes 
have then been mosaiced into a master datacube cov- 
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ering the entire ring and its surroundings. 

For the stellar observations, columns containing 
high signal-to-noise spectra are extracted and summed; 
Gaussian profiles are fitted to the hydrogen recombi- 
nation lines (e.g., Pa/3 at 1.28 ^m, Br7 at 2.16 /^m) 
then used to fill in this absorption; and finally they 
are divided by the Planck function appropriate to a 
blackbody of the same temperature as the star. This 
yields response spectra which are divided into each 
spectrum of the J-band and K-hand master datacubes 
to compensate for variations in atmospheric transmis- 
sion with wavelength. 

Inspection of the datacubes does indicate the pres- 
ence of extended line emission at the redshifted wave- 
lengths of Pa/3, Br7, and H2 1-0 S(l). However, the 
net signal-to-noise in the if-band master datacube is 
significantly poorer than the J-band, and inadequate 
for our purposes. We therefore have to rely solely on 
the J-band datacube for our analysis. The equiva- 
lent width of the Pa/3 emission line has been measured 
in the datacube by integrating the flux within the line 
and dividing by the interpolated continuum level. Fig- 
ure [S] shows a map of the Pa/3 equivalent width after 
binning x3 to match the pixel scale of IRIS 2, overlaid 
with contours of the (J — H) colour from Figure [J] 
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20''03°29.40' 29.20= 29.00" 28.80= 28.60= 
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Figure 6: Contours of (J — H) colour in the nu- 
cleus of IC 4933 from IRIS2, superimposed on a 
map of the Pa/3 equivalent width from GNIRS 
IFU spectroscopy. The contours correspond to 
(J - Jf) = (0.71,0.73,0.75,0.77,0.79), while the 
equivalent width in units of A is indicated by the 
scale at right. 



By contrast with the broad-band colour maps, the 
Pa/3 map shows a clear bifold symmetry, with the equiv- 
alent width reaching a maximum of 14 ± 2 A almost 
due North and South around the ring, and a minimum 



of 5.5 ± 0.5 A in the East and West. Figure [7] shows 
the evolution of Pa/3 equivalent width with time for 
the IB and CSF Starburst99 models of Section ETT] 

The IB model reaches the Pa/3 levels of the North/South 
quadrants 6.6 Myr after the onset of the burst, and 
those of the East/West quadrants a further 0.3 Myr 
later. This is perfectly consistent with the IB colour 
evolution described by Fig. |4l with azimuthal varia- 
tions in extinction accounting for the morphological 
differences between the nuclear ring as seen in (J — H) 
and in Pa/3 (Fig. |6]). Indeed, the extinctions inferred 
for each quadrant from their derived ages and location 
in Fig. |4] are found to range from Ay ~ 2 mag in the 
E quadrant, up to 3.5 mag in the N quadrant. 

Although this plot covers only the first 10 Myr, 
even 1 Gyr after the onset of CSF the Pa/3 equiva- 
lent width is still 50% higher than the peak values ob- 
served anywhere around the ring. While it is conceiv- 
able that the observed Pa/3 equivalent widths are un- 
derestimated due to inadequate allowance for an older 
background stellar population contaminating the con- 
tinuum fluxes, the very gradual decline in Pa/3 equiv- 
alent width over time in the CSF model would make 
the East/West quadrants at least 0.5 Gyr older than 
the North/South quadrants. We can conceive of no 
scenario under which two populations of such differing 
ages could remain spatially distinct for so long in such 
a compact, but kinematically active region and there- 
fore discount purely continuous star formation in the 
circumnuclear ring of IC 4933 as playing any signifi- 
cant role. 
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Figure 7: The evolution with time of the Pa/3 
equivalent width for the same Starburst99 IB 
and CSF models as in Fig. The upper and 
lower shaded regions mark the observed range of 
equivalent widths within the North/South, and 
East/ West quadrants respectively, of the IC 4933 
nuclear ring. 



www.pu blish . csiro. au /journals/ pasa 



7 



3 Discussion 

The IB and CSF Starburst99 models represent the 
two extremes of star formation, and as stellar popu- 
lat ion model fi t s to o ptical line indices of nuclear rings 
bv lSarzi et alj (|2007l ) has indicated, the truth is more 
likely to be a series of multiple, episodic bursts super- 
imposed on an old underlying bulge population. Nev- 
ertheless, the age differences around the ring are un- 
mistakable. IC 4933 thus joins a small, but growing 
list of galaxies whose nuclear rings exhibit evidence of 
a clear bipo lar azimuthal age gradient such as MlOO 
(|Rvder et a l. 2001; AUard ct al. 2006), NGC 1343, and 
NGC 1530 (Mazzuca ct al. 2008). Many more galax- 
ies show some fo rm of azimuthal age gradient at least , 
like that in M83 (iPuxlev et al.ll 19971 : iHarris et all200ll : 



Knapen et al.ll2010l '). including 60% of the sample of 
Boker et alj l|2008t l observed in the NIR with the SIN- 



FONI IFU on the V LT, and 40% of the sample of 
Im azzuca et all (|2008l) ob served in Hq. 

iMazzuca et al.l l|2008l ) found no clear relationship 
between the presence or absence of an age gradient, 
and the morphology of either the ring or the host 
galaxy. They did note that the rings in their sam- 
ple that show no clear age gradient are also among 
the oldest (>10 Myr) star-forming rings. Thus it may 
be that with the passage of time and several complete 
orbits of the aging clusters around the ring that a well- 
organised, sequential star formation pat tern may sim- 
ply be smoothed over (|Boker et al.ll2008l 'l. We do know 
that the starburst in IC 4933 is comparatively young 
and of short duration, which may partly account for 
why it still has such a well-defined bipolar age gra- 
dient. A weak or intermittent gas inflow rate from 
the bar into the ring may also yield a more stochastic 
age distribution of clusters around the ring, particu- 
larly if star formatio n only occurs above some criti- 
cal threshold density ([Elmegreenlll997l : iMazzuca et al.l 
[2008; Boker et al. 20081'). 

IMazzuca et al.l (|2008l l plotted the angular offset 
between the bar major axis position angle and the 
azimuth of the youngest Hll region around each half 
of their nuclear rings. They found a broad distri- 
bution with an apparent peak around 90°, suggest- 
ing that soon after inflowing gas reaches the "contact 
points" between the offset dust lanes and the nuclear 
ring (commonly found perpendicular to the bar major 
ax;is) it is converted to stars. We have used the ellipse 
surface photometry task within iraf to derive the bar 
position angle. Over radial distances between 10 and 
18 arcsec where the bar dominates, the position angle 
is (32 ±2)°. The two Pa/3 peaks in Figure |S] have po- 
sition angles of 182° and 350°, so are displaced some 
30° — 40° from the bar position angle. The spatial 
resolution of Figure [S] is insufficient to discern if the 
contact points in IC 4933 are indeed at 90° to the bar 
major axis, or could in fact be much closer to these 
youngest Hll regions. 

The GNIRS Pa/3 datacube allows us to extract a 
crude velocity field which yields a deprojected rota- 
tion velocity around the ring of 220 kms~^ with the 
northeast side of the disk approaching us. Assuming 
a bar pattern speed in the region of 25 kms^^ kpc~^ 



(jRvder et al.lll996l l. newly-formed clusters would then 
drift away from the contact point at ~200 kms"'^. At 
that velocity, one complete orbit of the ring takes just 
23 Myr. So we should perhaps not be too surprised 
to find relatively young clusters so far from the in- 
ferred current location of the contact points, where 
star formation could be happening right now but will 
not become apparent even in the NIR for a few Myr. 

Al ternatively, as shown by the simulations of lHeller fc Shlosmanl 
(|2006l ). an oval (rather than circular) nuclear ring can 
result in contact points which are not perpendicular to 
the bar major axis. It is notable that the nuclear ring 
seen in the colour maps and in Pa/3 is quite circular, 
whereas IC 4933 appears to be inclined to the line-of- 
sight by 37° ( )Lauberts fc Valenti in 1989). Thus, fore- 
shortening of an eccentric nuclear ring could result in 
it appearing circular in the plane of the sky, and it may 
in fact be that the contact points of the dust lanes with 
the ring do coincide with the youngest H ll regions. 

Interestingly, the blue 'hot spots' in the (it' — K) 
image (Figure O also show large spatial offsets from 
the Pa/3 equivalent width peaks. The Starburst99 IB 
model in Figures|3]and[7]predicts (in the absence of any 
reddening) that (U — K) becomes bluer by 1 mag in the 
first 4 — 5 Myr, then reddens rapidly by 4 mag in the 
next 3 Myr and stays within 1 mag of this thereafter. 
For these hot spots to appear as blue as they do, in the 
presence of dust, would require that they be no older 
than 5 Myr. This is inconsistent with the chronology 
derived from Pa/3 equivalent widths, which are far less 
susceptible to reddening, and again highlights the risk 
in relying solely on broadband colour measurements for 
deriving star formation histories on short (<10 Myr) 
timescales. 



4 Conclusions 

Motivated by the fact that barred spiral galaxies known 
to possess large-scale inner rings are more likely to host 
nuclear rings as well, we obtained J, H, and Ks imag- 
ing of the nearby SB(rs)bc galaxy IC 4933 with IRIS2 
on the AAT, and it' imaging with CMOS on Gemini 
South. Although the reduced images show only the 
bulge, bar, and arms, an almost circular ring 5 arc- 
sec (1.5 kpc) in diameter is seen in the colour index 
maps. Comparison of the ring NIR colours with evo- 
lutionary tracks from Starburst99 models is unable to 
distinguish between a burst population just 5 Myr old, 
and less-reddened continuous star formation over the 
past 10-100 Myr. 

The equivalent width of the Pa/3 line is much less 
susceptible to reddening effects. We have mosaiced 
four pointings of the GNIRS Integral Field Unit in or- 
der to map the Pa/3 equivalent width around the ring. 
We find clear peaks in the equivalent width in the 
northern and southern sectors which are most easily 
accounted for by short-lived bursts of star formation 
6.6 Myr ago rather than by continuous star formation, 
with the eastern and western sectors appearing to be 
0.3 Myr older. IC 4933 is the fourth galaxy seen to 
have evidence of such a well-defined bipolar azimuthal 
age gradient. Although the actual location of the con- 
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tact points of the offset dust lanes in the bar with the 
nuclear ring cannot be clearly discerned in our images, 
the significant azimuthal offset between the youngest 
Hll regions in the ring and the bar major/minor ajces 
are consistent with this nuclear ring being intrinsically 
oval in shape but fortuitously appearing circular from 
our perspective. Ongoing and future NIR imaging sur- 
veys such as those conducted with UKIRT and VISTA 
can be expected to yield many more such examples of 
nuclear rings, but only a comprehensive NIR integral 
field spectroscopic followup can unambiguously define 
the star formation history around these rings. 
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